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ABSTRACT 

In quintessence models, the dark energy content of the universe is described by 
a slowly rolling scalar field whose pressure and energy density obey an equa- 
tion of state of the form p = wp; w is in general a function of time such that 
w < —1/3, in order to drive the observed acceleration of the Universe today. 
The cosmological constant model (ACDM) corresponds to the limiting case 
w = —1. In this paper, we explore the prospects of using the Lyman-a for- 
est to constrain w, using semi-analytical techniques to model the intergalactic 
medium (IGM). A different value of w changes both the growth factor and 
the Hubble parameter as a function of time. The resulting change in the op- 
tical depth distribution affects the optical depth power spectrum, the number 
of regions of high transmission per unit redshift and the cross-correlation co- 
efhcient of spectra of quasar pairs. These can be detected in current data, 
provided we have independent estimates of the thermal state of the IGM, its 
ionization parameter and the baryon density. 

Key words: Cosmology: theory - intergalactic medium - large-scale structure 
of universe - quasars: absorption lines 



1 INTRODUCTION 

Observed cosmic microwave background anisotropics 
convincingly demonstrate that the universe is spatially 
flat (de Bernardis et al. 2002; Netterfield et al. 2002; 
Pryke et al. 2002). The luminosity-distance relation, 
as determined from high redshift type la Supernovae 
(Garnavich et al. 1998; Riess et al. 1998; Perlmutter et 
al. 1999), requires a spatially flat universe to be cur- 
rently dominated by some type of dark energy, which 
is nearly homogeneous and has negative pressure, such 
as for example a cosmological constant. Several inde- 
pendent lines of argument also favour a low-density, 
vacuum-energy dominated universe, for example the 
abundance of clusters (Bahcall et al. 2003) and their 
X-ray properties (Allen et al. 2002), the clustering of 
galaxies (Efstathiou et al. 2002; Verde et al. 2002), es- 
timates of the age of the universe and the current value 
of the Hubble parameter. 

The energy density associated with the cosmo- 
logical constant may actually decay in time. Such 
quintessence models (e.g. Caldwell et al. 1998; hence- 



forth QCDM models) have an equation of state p = w p, 
with w < —1/3, where w — —1 corresponds to the more 
familiar cosmological constant model. Recently, a num- 
ber of observational tests have been proposed to mea- 
sure w and its redshift dependence (Kujat et al. 2001; 
Wang & Garnavich 2001; Matsubara & Szalay 2002; 
Gerke & Efstathiou 2002). 

In this Letter we focus on the prospects of using 
the Lyman-a forest to constrain w. The Lyman-a for- 
est region in quasar (QSO) spectra contains hundreds 
of hydrogen Lyman-a absorption lines (see Ranch 1998 
for a review) , most of which are produced by small den- 
sity fluctuations in the intervening intergalactic medium 
(Gen et al. 1994; Bi & Davidsen 1997, hereafter BD97). 
Since these structures are still reasonably linear, they 
are good probes of the underlying large-scale matter 
distribution, and hence are sensitive to the cosmologi- 
cal model. For example, Hui et al. (1999) and McDon- 
ald (2001) suggested to apply the Alcock-Paczinski test 
(Alcock & Paczinski 1979) to estimate the cosmological 
constant and Croft et al. (2002) constrained the dark 
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matter power spectrum on large scales. The Lyman- 
ot forest is promising in that it can be used over a larger 
redshift range than most other potential probes of w. 

Here we use semi-analytical models of the Lyman- 
a forest, which have boon shown to be successful in re- 
producing reasonably well most of its observed proper- 
ties, such as the column density distribution function, 
the line-widths distribution and the number of lines 
per unit redshift (BD97; Viel et al. 2002a, 2002b; Roy 
Choudury et al, 2001). Hydrodynarnical simulations arc 
required for more accurate predictions on smaller scales 
(e.g. Theuns et al. 1998, 2002; Dave et al. 1999; Bryan et 
al. 1999). Wc examine several statistics of the flux spec- 
trum, such as the probability distribution function of 
the optical depth, the optical depth power spectrum, the 
number of underdense regions per unit redshift and the 
cross-correlation coefficient of spectra of quasar pairs. 
The first two statistics arc not directly observable but 
they can be recovered from the flux distribution by look- 
ing at higher order lines through pixel optical depths 
techniques (Sect. 3). Section 2 briefly reviews the prop- 
erties of quintessence models and describes our simu- 
lation techniques. Section 3 contains our results and a 
brief discussion. 



2 SIMULATIONS OF Lyman-Q SPECTRA 
IN QUINTESSENCE MODELS 

The quintessence parameter w influences the shape 
of the linear matter power spectrum P{k, z) oc 
k'^ T^{k, z) D'^{z), the linear growth factor of density 
perturbations D{z) and the Hubble parameter H{z). 
Here k is the wavenumber, T the transfer function 
and n the spectral index. Useful expressions for the 
growth factor of density perturbations and of linear pe- 
culiar velocities can be found in Lahav et al. (1991), 
for ACDM, and in Wang & Steinhardt (1998), for 
QCDM. The Hubble parameter at redshift z is H{z) = 

Ho [n^(i + zf + nfc(i + zf + noEii + zf+^^] 

with Ho = 100 /i km s~^ the present-day Hubble pa- 
rameter and Q,k — 1 — fim — Odb. Odb represents the 
energy density of the cosmological constant {Q.a) or the 
quintessence (Oq), in units of the critical density. 

In Figure 1 we plot T{k,z), D(z) and the Hub- 
ble parameter H{z), for several values of w, here as- 
sumed to be constant in time, for simplicity, using the 
fits of Ma et al. (1999). If we compare the transfer func- 
tion of QCDM and ACDM at « = 2 one can see that 
there are differences only on very large scales, k ~ 0.01 h 
Mpc~^. The Hubble parameter in QCDM models with 
w = -0.8, -0.6 and -0.4 differs by 6, 14 and 25 per 
cent at redshifts z ~ 0.75, « ~ 1 and z ~ 1.2 from those 
of the ACDM model (panel b). For those values of w, 
the growth factors D{z) (panel c) at z ~ 2 are larger 
than for the ACDM model by 5, 15 and 30 per cent 
respectively. 

The semi-analytical model we use to generate 
Lyman-a forest spectra is based on the approxima- 
tions introduced by Bi and collaborators (BD97; see 
Viel et al. (2002a) for more details). The basic as- 
sumption is that the low-column density Lya forest is 



produced by smooth fluctuations in the intergalactic 
medium (ICM) , which arise as a result of gravitational 
growth of perturbations (e.g. Schaye 2001). Briefly, we 
start by generating correlated Gaussian random fields to 
represent the linear density and peculiar velocity dark 
matter fields along a sight-line, for a given linear mat- 
ter power spectrum. The linear density perturbations 
in the IGM So^^{'x.,z) are related to those in the un- 
derlying dark matter by a convolution, Sl'^^ (k, z) = 
WiGM{k, z)5^'^'^{'k, z). For the smoothing kernel W we 
use the Gaussian filter Wigm = exp(— fe/fe/)^ (Gnedin 
et al. 2002), which is a good approximation to model the 
effects of gas pressure in the linear regime. The smooth- 
ing scale kf 18 related to the Jeans length in a way which 
depends on the thermal history of the IGM (Gnedin & 
Hui 1998; Matarrese & Mohayaee 2002). 

To account for the fact that the IGM which pro- 
duces most of the Lyman-a forest is actually mildly non- 
linear, we follow BD97 and adopt a simple lognormal 
model (Coles & Jones 1991) for the IGM local density. 
This is a strong assumption but we think that at these 
redshift and for the large scale properties investigated 
here the log-normal model can be a good approxima- 
tion to the more accurate hydro-dynamical simulations. 
We assume the gas to be in photoionization equilib- 
rium with am imposed uniform UV-background, which 
we scale so that the mock spectra have the same effec- 
tive optical depth r^a as observed (see Kim et al. 2002 
and Bernardi et al. 2002, for recent determinations of 
Teff{z)). We use the values of Schaye et al. (2000) for 
the normalization To and slope 7 of the temperature- 
density relation, T = To(p/(/9))^"\ in the IGM. Fi- 
nally, we artificially broaden the spectra to a resolu- 
tion of FWHM = 6.6 km s"^ to mimic high quality 
VLT/UVES or Keck/HIRES spectra, and add Gaus- 
sian noise with signal to noise of 50 (sec e.g. Theuns et 
al. 2000 for more details). Spectra generated with this 
procedure have been shown to produce flux probabil- 
ity and line-width distribution functions in reasonable 
agreement with observations (BD97). 

We will compare the Lyman-Q forest in a ACDM 
cosmological model with three different QCDM models, 
with w = —0.4, —0.6, —0.8, respectively. For all of these 
models, we assume that the cosmological parameters are 
as = 0.7, Qm = 0.3, ft = 0.67, fi,, ft^ = 0.020, Qde = 
0.7. We simulate 10 different randomly generated spec- 
tra, with the same set of random phases for each spec- 
trum, in the redshift range 1.8 < z < 2.2 for each model, 
assuming r^ff = 0.15 (Kim et al. 2002), 7 — 1.3 and 
To = lO"'^ K. The fihering scale is kf ~ 1.7 fcj (with 
kj the Jeans wave-number), which is a reasonable es- 
timate if HI reionization takes place at a ~ 7 (Gnedin 
& Hui 1998). The spectra are approximately 40, 000 km 
s~^ long. In addition, we simulate spectra of QSO pairs 
with a given angular separation in the range 10-90 arc- 
sec, using the procedure described in Viel et al. (2002a). 
The cross-correlation coefficient spectra of pairs is esti- 
mated using the definition of Viel et al. (2002a) from 
a set of 8 pairs for each separation. In panel d) of Fig- 
ure 1 we show, for a qualitative comparison, a chunk of 
two simulated spectra for the ACDM and QCDM, with 
w = —0.4, model. We have decided to focus our analysis 
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Figure 1. Panel a): ratio of transfer functions Tqcdm /Tacdm ctt z=2. Panel b): ratio of Hubble parameters Hqcdm /H^cdm 
as a function of redsliift. Panel c): ratio of growth factors Dqqdm / D\cDM a function of redshift. The three curves shown 
are for three quintessence models with w = —0.4 (continuous line), w = —0.6 (dotted line) and w = —0.8 (dashed line). Panel 
d): example of simulated spectra, it is shown ACDM spectra (dot-dashed line) and the QCDM with w = —0.4 (continuous line). 



around z ~ 2 in view of a comparison with observations 
which will be made in a future work. 



3 RESULTS 

In this section we will study the effects of dark en- 
ergy on the simulated spectra. The differences between 
the QCDM and ACDM model, which could be tested 
through Lyman-ci absorptions, are the following: the lin- 
ear dark matter power spectrum P{k,z); the evolution 
of the cosmological parameters and in particular of the 
Hubble parameter H{z); the growth factors D{z) of the 
linear density perturbations. 

The simulated Lyman-a flux power spectrum 
agrees reasonably well with the linear dark matter power 
spectrum at large scales, while on smaller scales non- 
linear effects, thermal broadening and noise can produce 
differences. This is roughly true even if we do not cor- 
rect for redshift-space distortions and compare directly 
the flux power spectrum of simulated absorption spec- 
tra and the linear dark matter power spectrum. From 
Figure 1 (panel a) we can see that the QCDM trans- 
fer function differs from the ACDM one only on scales 
larger than 100 comoving Mpc. Such scales are compara- 
ble to the extent of the Lyman-a forest region in a QSO 
spectrum, and hence this signature of w —1 would be 
difficult to detect since it would require a large number 
of spectra to estimate the power on these scales. More 
importantly, in high resolution spectra every attempt to 
recover the power spectrum on scales larger than 10 co- 
moving Mpc is challenging due to the continuum fitting 
uncertainties (Croft et al. 1999). 

A second difference between QCDM and ACDM is 
the different evolution of the cosmological parameters, 
such as H{z). We have explored the effect of a different 
Hubble parameter and different growth factors values 
at redshift z = 2 on the simulated Lyman-o? optical 
depth, by fixing the remaining parameters. The effect 



is clear: a larger Hubble parameter produces in gen- 
eral a broader optical depth distribution with the peak 
shifted to lower values. In fact, with a higher value of 
the Hubble parameter, the same real-space absorber in- 
tersected by the sight-line will affect a larger region of 
the redshift-space spectrum. 

In the same way, we explore the effect of differ- 
ent growth factors. If the growth factor is increased, the 
(non-linear) log-normal mapping acts in such a way that 
typically the low-density regions become less dense and 
the high-density regions more dense. Thus, with a larger 
growth factor, as the one predicted for QCDM models 
compared to ACDM ones (Figure 1, panel c), the simu- 
lated optical depth shifts to smaller values and its proba- 
bility distribution function (pdf) becomes broader. This 
effect can be seen also in the simulated fiux in panel d) 
of Figure 1, for example at ~ 6000 km s~^ this is par- 
ticularly evident. 

Our first conclusion is that a larger growth factor 
produces the same effect of a higher value of the Hub- 
ble parameter: a shift of the simulated Lyman-a optical 
depth to smaller values and a broader distribution, if 
the other parameters remain fixed. 

To see these effects combined together, we plot in 
panel a) of Figure 2 the probability distribution func- 
tion of the simulated optical depth for ACDM (dot- 
dashed line), QCDM with w — —0.4 (continuous line), 
w — —0.6 (dotted line) and w — —0.8 (dashed line). 
We choose to plot the probability distribution function 
also for T values smaller than those which can actu- 
ally be recovered by standard techniques to better ap- 
preciate the differences between cosmological models. 
We can see that the simulated optical depth pdf re- 
sults in different distributions for the four cosmolog- 
ical models. Since the flux is the observed quantity, 
the inferred optical depth can be trusted only in the 
range —1.6;$ log(T)^ 0.5 (corresponding to fluxes in 
the range 0.04-0.97). With a more sophisticated anal- 
ysis, based on pixel optical depth techniques, this range 
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Figure 2. Panel a); probability distribution function of the optical depth (panel a) for QCDM model with ui = —0.4 (continuous 
line), UI = —0.6 (dotted line), w = —0.8 (dashed line) and ACDM model (dot-dashed line). The interval in which the optical 
depth can be recovered directly from the flux is represented by the double arrow (a larger range can be used with pixel optical 
depth technique, see text). Panel b): power spectrum of the optical depth; error bars represent the error of the mean. 



can be considerably extended (Aguirre et al. 2002 and 
references therein) . In panel b) of Figure 2 we plot the 
optical depth power spectrum for the QCDM and the 
ACDM model, with Poissonian error bars. We can see 
that on scales logfc(kms~^) < —1.5 the four cosmolog- 
ical models are distinguishable: while the shape is very 
similar, the normalization is different even if we choose 
the same value for r^ff. This is the effect of different 
amounts of non linearity in the four cosmological mod- 
els, which influences the normalization of the curves. 

In Figure 3 we plot two more observational quan- 
tities: the number of underdense regions per unit red- 
shift (panel a) and the cross-correlation coefficient x es- 
timated from the spectra of QSO pairs at different sep- 
arations (panel b). The definition of underdense region 
adopted here is the following: if we are searching for an 
underdense region of size Aw, we preliminarily smooth 
the simulated spectra on this scale (top-hat smoothing) 
and then we define as an underdensity a region for which 
the flux is larger than the mean flux at that redshift, es- 
timated through Te//. One can see that the predicted 
number of underdense regions per unit redshift is dif- 
ferent for the four cosmological models. The number of 
underdense regions of intermediate sizes, between 700 
km s^^ and 2700 km s~^, which correspond to sizes in 
the range 10-40 comoving Mpc in a ACDM universe, 
is larger in the ACDM model than in the QCDM ones. 
This is due to the combined effects of a different Hubble 
parameters and a different growth factor for density per- 
turbations, as we saw in panel a) of Figure 2, which pro- 
duce more rarefled underdense region and denser high 
density regions in QCDM models than in ACDM mod- 
els. In particular, we have found that the number of 
underdense regions of size 2500 km s~^ for ACDM can 
be a factor of 2 and 1.5 larger than in QCDM model 
with w — —0.4 and w = —0.8 at z = 2, respectively. 

The cross-correlation coefficient x. is plotted in 
panel b) of Figure 3, this number quantifies the amount 
of 'coherence' in the transverse direction between spec- 
tra of QSO pairs. The information contained in the 



transverse direction could be powerful and has already 
been analysed in the past (Miralda-Escude et al. 1996; 
Viel et al. 2002a). This test is in a sense similar to 
the more sophisticated Alcock-Paczinski test (Hui et 
al. 1999; McDonald & Miralda-Escude 1999; McDon- 
ald 2001) and aims at discriminating between different 
cosmological models. Here, we have simulated 8 QSO 
pairs at 2 2 for five separations (10, 30, 50, 70 and 
90 arcsec) and computed the mean value and the error 
of the mean for this coefficient at each separation. The 
four models are very similar and the error bars overlap. 
Only at angular scales ^ 20 arcsec there are some differ- 
ences between the models. However, the number of pairs 
considered here is marginally sufficient to discriminate 
between a ACDM model and a QCDM with w = -0.4. 
We have estimated that 25 pairs with these small sep- 
arations are needed to distinguish between the ACDM 
model and a QCDM with w = —0.6 at a 3cr level. 

Dark energy in the form of quintessence in the 
Lyman-a forest seems to be not easily distinguishable 
from a cosmological constant. Even if the thermal state 
of the IGM is completely known, the uncertainties on 
the values of the baryon fraction (Buries et al. 2001) 
and especially on the ionization parameter (Scott et 
al. 2000) , are so large that it would be difficult to recover 
the evolution of H{z). We have found that a QCDM 
model with w = —0.4 requires an ionization parameter 
a factor of 3 lower than in a ACDM model (while for 
w — —0.8 it is a factor of 1.25 lower), in order to match 
observed r^f / . Uncertainties on the value of r^f / are of 
the order of 40%, but in the near future, with larger 
data set available, this uncertainty will become signifi- 
cantly smaller. In principle, an accurate determination 
of Teff, with a precision of 10 — 20% could already be 
used to discriminate between QCDM with w — —0.4 
and ACDM models, since the present uncertainties on 
the ionization parameter are of the order of 50%. Statis- 
tics based on the optical depth like the pdf or the power 
spectrum, instead of the flux, can be much more useful 
in discriminating between different cosmological mod- 
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Figure 3. Number of underdense regions per unit redshift as a function of void size in km s^^ (panel a) for the three QCDM 
models with w = —0.4 (continuous line), ui = —0.6 (dotted line), w = —0.8 (dashed line) and ACDM model (dot-dashed 
line). The simulated sample consists of 10 spectra, spanning the redshift range 1.8 < z < 2.2 for each cosmological model. 
Cross-correlation coefficient, x, as a function of QSO pair separation in arcsec (panel b). The simulated sample consists of 8 
QSO pairs for each separation shown. Poissonian error bars in panel a), shown only for two different models for clarity. In panel 
b) the error bars are the errors of the mean value. 



els. This is due to the fact that the optical depth is very 
sensitive to the amount of non-linearity introduced in 
the different cosmological models. However, in order to 
do that, it will be necessary to use pixel optical depth 
techniques, using higher order Lyman lines, to get reli- 
able estimates of the optical depths both in saturated 
and in high transmissivity regions. 
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In summary, we have presented a preliminary anal- 
ysis of simulated spectra in QCDM models and in a 
ACDM model based on semi-analytical models. The 
main differences can be found in the different evolution 
oi H{z) and of the linear growth factor of density pertur- 
bations D{z). Both these effects act in a similar way. If 
we fix the remaining parameters and compare the prob- 
ability distribution function of Lyman-a optical depth 
between QCDM and ACDM we have found that the 
distribution in the QCDM model is broader than in the 
ACDM case. This means that in general high transmis- 
sion regions are more rarefied and high density regions 
are denser in QCDM than in a ACDM universe. The op- 
tical depth power spectrum is different on intermediate 
and large scales, even if the normalization Te// is the 
same, because of the different amounts of non linearity 
predicted for the models. Detecting the underdense re- 
gions could be promising. In fact, we have found that the 
number of intermediate size underdense regions, with 
sizes of the order of 10 comoving Mpc, is larger in a 
ACDM model than in QCDM ones. Another statistics 
we have proposed concerns the the cross-correlation co- 
efficient between the spectra of QSO pairs. However, in 
this case, the differences can be appreciated only for 
separations smaller than 20 arcsec and the number of 
pairs should be larger than 20 to distinguish between 
QCDM and ACDM cosmological models at a Scr level. 
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